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We reinvestigate the competition between the Mott and the Anderson insulator state in a one- 
dimensional disordered fermionic system by a combination of instanton and renormalization group 
methods. Tracing back both the compressibility and the ac-conductivity to a vanishing kink energy 
of the electronic displacement field we do not find any indication for the existence of an intermediate 
(Mott glass) phase. 
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Introduction. Two non-trivial routes to insulating be- 
havior in solids are connected with the names of Ander- 
son and Mott. Non-interacting electrons in disordered 
solids are always localized in one and two space dimen- 
sions resulting from a combination of classical localiza- 
tion and quantum- mechanical interference [J Q. This 
Anderson insulator (AI) is characterized by both a finite 
ac conductivity at low frequencies and a finite compress- 
ibility Interaction between electrons reduces the ef- 
fect of disorder and may lead to metallic behavior in two 
dimensions In Mott insulators, on the other hand, 
insulating behavior results from the blocking of sites by 
repulsive interaction between electrons [5| and hence are 
dominated by correlation effects. The Mott-insulating 
(MI) phase is incompressible and has a finite gap in the 
conductivity. A natural question is then: what happens 
in systems when both disorder and interactions are non- 
negligible? Is there, depending on the strength of in- 
teraction or disorder, a single transition between these 
two phases or is the scenario more complex? A particu- 
larly interesting case to answer these questions is that of 
electrons in one space dimension: interactions are very 
strong, destroying Fermi liquid behavior, but interac- 
tions can be partially incorporated into the harmonic 
(bosonized) theory [6J]. Additional motivation to study 
this case comes from its physical realization in quantum 
wires Q and ultra-cold gases [||. 

One-dimensional systems of this type have indeed been 
considered in a number of publications. Early work by 
Ma using real-space renormalization group (RG) for 
the disordered ID Hubbard model suggests a direct tran- 
sition from an AI to an MI phase. This result contrasts 
with more recent work (see e.g. 13; U| and references 
therein) where a new type of order, different from the MI 
and AI phase, was found. In particular in [llj the exis- 
tence of a new Mott glass (MG) phase was postulated 
which is supposed to be incompressible but has no gap 
in the optical conductivity. 

In this paper we reinvestigate this problem by relating 
both the compressibility and the low frequency conduc- 
tivity to the energy of kinks in the displacement pattern 
of the bosonized electrons, resorting to ideas known from 
the discussion of the roughening transition Using 



an instanton calculation of the conductivity 1JJ, 14 1 we 
show that an incompressible phase has always a gap in 
the low frequency conductivity, excluding the possibil- 
ity of a MG phase. A number of further consideration 
support this finding. 

The model. Following Haldane [15| we relate the mass 
density p(x) of the electrons to their displacement field 
(p(x): p(x) = po + ±d x ip + p cos(2ip + 2k F x) + . . .. 
Throughout the paper we will use e = h = 1. The 
ground state energy of the system follows from E = 
— liniT-^o TTn J D(pe~ s where the action is @, [ll[ 
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Here we introduced dimensionless space and imaginary 
time coordinates by the transformation Ax — ► x and 
Avt — > t where v is the plasmon velocity and A a 
large momentum cut-off. At = v/T and L and de- 
note the thermal de Broglie wavelength of the plas- 
mons and the system size, respectively, w denotes the 
strength of the Umklapp scattering term (or the strength 
of a commensurate potential) and F the external driv- 
ing force. fi(x) and £(x) result from the coupling of 
the random impurity potential to the long wave length 
and the periodic part of the density, respectively, with 
(C(x)C*(x')) = u 2 S(x-x'), (p{x)p{x')) =aS(x-x'), all 
other correlators vanish. Finally we add a dissipative ac- 
tion Sdiss describing an Ohmic resistance which may e.g. 
result from the interaction with an external gate [161 ] . 

Rigidities. We begin with a discussion of the gener- 
alized rigidities, which are related to the compressibility 
and the conductivity of the system. We first consider the 
application of a fixed strain by imposing the bound- 
ary conditions ip(0,r) — and ip(L,r) — irdLA. The 
boundary condition in the r-direction is assumed to be 
periodic. For $ <C 1 and L — * oo the corresponding in- 
crease AEoitf, 0) = Eoi'd, 0)-E (0, 0) of the ground state 
energy Eq^, 0) is clearly an even but not necessarily an- 
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FIG. 1: Left: Kink and anti-kink in the displacement profile 
tp(x). The thin lines represent the minima of the potential 
energy in the absence of the driving force, the bold line one 
metastable state and the dashed line the instanton configu- 
ration, respectively. Applying a fixed strain to the system 
only kinks (or anti-kinks) are enforced in the system. Right: 
schematic phase diagram with w c ~ a 2 . 
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The r.h.s. of this relation has to be understood as follows: 
if £ x =0 then the stiffness T x = A 2 /k describes the 
response to the twisted boundary conditions, T" 1 is the 
isothermal compressibility [l3(. In this case the change 
of tp is spread over the whole sample. If however T x 
diverges, i.e. the system becomes incompressible, then 
the kink energy T, x /A is non-zero. In this case the change 
of if from to it occurs in a narrow kink region of width 
£ much smaller than L. Creating a kink corresponds to 
adding (or removing) an electron at the kink position. A 
non-zero kink energy resembles the step free energy of a 
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surface below the roughening transition Yt 

In a similar manner we can apply non-trivial boundary 
conditions in the r-direction by choosing (p(x, 0) = and 
tp(x,L T ) = 7tjL T /v. This corresponds to imposing an 
external current j = (d T <p) /tt at x = and x = L: 
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if E r = 0, (4) 
if T- 1 = 0. (5) 
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1/T T = kv 2 . D determines the Drude peak of the con- 
ductivity a(uj) = D5(uj) + cr rog (u;) @, In Lorentz 
invariant systems (like the MI) T x — (Av) 2 T T and 
£ x = AwE T , provided they are finite. 

So far we assumed that T x / T , £ x / r are self-averaging 
if L — > oo. In the same way we may introduce local 
rigidities by applying twisted boundary conditions over 
a large but finite interval [x, x + L X A], A -1 « 4 < L. 
The result for T x / T and E x / r will then depend in general 
on the size L x of the interval , i.e. E x / T — > Ti x / T (L X ) . 

Conductivity. In cases where D vanishes, the field or 
frequency dependent conductivity a(F 1 uj) may still be 
non-zero, provided Forw are finite. The non-linear dc- 
conductivity a(F, 0) follows from an instanton approach 
in which the current is calculated from the tunneling 
probability of electrons between consecutive metastable 



states (see Figfl]). The latter are the classical ground 
states in the absence of the driving force which follow 
from each other by a shift of a multiple of tt. This ap- 
proximate treatment is reliable for not too large K and 
small F. The conductivity lna(F, 0) = — 5i(L x ,s> £ T ,s) is 
then related to the action Si of a critical droplet (the in- 
stanton) of the new metastable state on the background 
of the old one [3, Hj- Here L XtS ,L TjS are the saddle 
point values of the extension of the instanton which fol- 
low from 

S*i« 2 ^ Lx) L T + L X (2 E r + ry In ^) - —FL X L T . (6) 
Av 4 irv 

For simplicity we assumed that the instanton is rectangu- 
lar, a choice which is motivated by the fact that the dis- 
order is correlated in time [l^ . E x /(Aw) and E T play the 
role of the surface tension of the instanton. The last term 
is the bulk contribution from the external field. The rj- 
term results from the dissipation which allows relaxation 
in the next metastable state [13j . 

The low frequency ac- conductivity a(0,u) results from 
the spontaneous tunneling processes between metastable 
states and their instanton configurations (1 and 2 in Fig. 
[TJ. Tunneling leads to a level splitting of the order 
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which has to match the energy ui of the external field [lj] . 
1/E T ~ K£ plays the role of the tunneling length, C is 
a numerical factor. This mechanism was first considered 
for non-interacting electrons by Mott [|| and extended 
to the interacting case via instantons in 14]. Thus, to 
get a non-zero <r(0, u>) for arbitrary low frequency u>, 
2T, X (L X ) < loA — > which requires 2Y, X {L X ) — > for 
a finite density of kink positions, which implies a finite 
compressibility. The distance L x between the sites in- 
volved in the tunneling is then L x (ui) » £K ln(l/(K^w)). 

Fixed points and phases. We come now to the discus- 
sion of the possible phases of model (JTJ) by attributing 
them to their RG fixed points (denoted by superscript 
*). Bare values will get a subscript o- For small u and w 
the lowest order RG-equations read [l9j |: 

^ = -K(au 2 + bw 2 ), ^ = a(l - cw 2 ) (8) 
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where / is the logarithm of the length scale and a, b, c are 
positive non-universal constants, z denotes the dynam- 
ical critical exponent which has to be determined from 
the fixed point condition. 
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Luttinger liquid (LL). The LL phase is characterized 



by u\ 



and hence E, 



0. 



1, 



K* L > 



and k* l = K l /{ttv* l ) > 0. The fixed point is reached for 
sufficiently large values of K and a. The long time and 
large scale behavior of the system is that of a clean LL 
characterized both by a finite compressibility k* l and a 
finite charge stiffness Dl = k* l v^. The dynamical con- 
ductivity is given by a reg — iD^j (ituS). The presence of 
the forward scattering term ~ \i{x)d x Lp does not change 
these results since it can be always removed by the trans- 
formation <p = (p — J* dx'[i(x'). 

Mott insulator (MI). Here K* M = k* m = u* M = a* M = 
but w* M ^> 1. Clearly the fixed point w* M is outside the 
applicability range of the ([8|)- (fTT|) but nevertheless some 
general properties of this phase can be concluded. The 
system is in the universality class of the 2-dimensional 
classical sine Gordon model which describes inter alia 
the MI to LL transition and the roughening transition 
of a 2-dimensional classical crystalline surface [13]. In 
the MI phase T x and T y diverge. The classical ground 
state is given by tp(x) — nix with n integer. The system 
is characterized by a finite kink energy E x ~ A(ko£m) _1 
where £m denotes the correlation length of the Mi-phase 
0. From © we get a(F,0) ~ exp(- F M / \K F)) where 
Fm ~ 47rS al E T /A ~ 1/( k o£m) is a characteristic depin- 
ning field 22j . According to (O the ac conductivity van- 
ishes for ijj < 2£a;/A. 

Anderson insulator (AI). Here w* A — K\ = but 
u* A 3> 1. k* l Rj kq is finite which is the result of the 
so-called statistical tilt symmetry [23[ corresponding to 
2 = 1. The fixed point Hamiltonian is in the univer- 
sality class of the 1 + 1 dimensional sine Gordon model 
with a random phase correlated in the r-direction. The 
transition to the LL phase occurs at K = K* A (u) > 3/2. 

Next we look at T, X (L X ) finite length scale £a <C L x <C 
L such that the parameters are close to their fixed point 
values. £a is the correlation length which diverges at the 
Beresinskii-Kosterlitz-Thouless (BKT) transition to the 
LL phase. To find the classical ground state of the sys- 
tem under periodic boundary conditions we have first to 
choose 2ifi + cxi — 2irni with m integer, and secondly 
the elastic term has to be minimized with respect to 
the ni. The subscript i refers to the sites of the lat- 
tice with spacing £4 and Q = |£de m \ The solution 
is ni = J2j<ii( a j ~ a i-i)/ 27I "]G EH- [x]g denotes the 
closest integer to x. The ground state is uniquely deter- 
mined by the aj apart from the pairs of sites (of mea- 



sure zero) at which aj 



±7r. At such pairs the 



ground state bifurcates since two solutions are possible. 
In the case of non- interacting electrons (K = 1) bifur- 
cation sites are states at the Fermi energy. For pairs 
at which aj — aj-i — ±ir + e, |e| < 1 we can go over 
to an excited state by creating a kink which costs at 
most the energy H x w cA/(k£a) iM- Those 'almost' 
bifurcating sites correspond to states close to the Fermi 
energy. The smallest e found with probability of order 



one in a sample of length L x is of the order £,/L x and 
hence S x « A/(L x kq), i.e. the kink energy vanishes for 
L x — > 00 [HI] and hence the system is compressible. 
Twisted boundary conditions in the y-direction give 



0) 



The non-linear conductivity is 
where Fa ~ 



given by o(F,uj = 0) - exp -y/F A /(FK 2 ) 

uM- This result can be understood as electron 
tunneling between sites at which Y, x ~ 1/L X . As ex- 
plained already a vanishing Y, x is also crucial for the 
existence of the low frequency conductivity a(Q,uj) ~ 
£m-F (wkqKZm ln(Ko^M^)) a s has been discussed in de- 
tail in [lj]. This result can be understood in terms of 
tunneling processes between rare positions at which the 
kink energies T, x (x) are much smaller than 1/(kqL x ). 



TABLE I: Properties of phases. 
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Mott- glass (MG). This new hypothetical phase was 
proposed in [111 ] to be characterized by a vanishing com- 
pressibility, Kq — 0, but a non-zero optical conductivity 
at low frequencies ll|. Since the phase is considered to 
be glassy, both fixed point values w G , ttg > 1. Similarly 
to the AI, the ground state can be found by minimiz- 
ing first the two backward scattering terms followed by 
minimization of the elastic energy. Although the ground 
state solution is now more involved than for the MI and 
AI case, for F = it is clearly periodic with period tt. 
As before kinks (or anti-kinks) with 5ip — ±ir allow the 
accommodation of twisted boundary conditions and the 
formation of instantons. A vanishing compressibility cor- 
responds to a finite kink energy T, x > C > which, 
according to leads to a gap in the ac- conductivity. 
Here we make the reasonable assumption that it is the 
instanton mechanism which dominates the low frequency 
response [20| . Thus, in a system with a non-zero cf(lu) 
for small uj also the compressibility has to be non-zero, 
contrary to the claims in [ll[ . 

Phase diagram. We come now to the discussion of the 
phase diagram of our model (JTJ) . From (110p follows that 
the random backward scattering term is generated by for- 
ward scattering and the commensurate pinning potential. 
Since a (I) = aoe 1 the two Eigenvalues Ai = 3 — 2K and 
A2 = 4 — 2K — 4a /n describing the RG-flow of u 2 and w 2 
around the the LL fixed point u* — w* — have oppo- 
site sign: u(l) increases whereas w(l) decreases. Thus the 
hypothetical MG phase, if it existed, cannot reach up to 



the point u = w = 0, in contrast to the findings in 11|. 
From this we conclude that for not too large values of Wq 
the AI phase is stable, as shown in Fig. [TJ To find the 
phase boundary to the MI phase we consider the stability 
of the MI phase with respect to the formation of a kink 
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by the disorder. To lowest order in the disorder we get 
for the kink energy in the MI phase 
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which gives for the phase boundary between the MI and 
the AI phase as depicted in Fig. [I] A similar result follows 
from the self-consistent harmonic approximation. 

So far we considered only typical disorder fluctua- 
tions. If \G(x)\ and n(x) are Gaussian distributed and 
rare events are taken into account, then Q12p remains 
valid with the replacements ctq — * &o ln(LA) and Uq —> 
Uq ln(LA). Thus the size of the MI phase is reduced but 
finite unless L — > oo. 

These findings are supported by the observation that 
the forward scattering term in (fT]) can be removed 
from the action by the transformation ip(x) = ip(x) — 
dx' ' n{x') j24j. On sufficiently large scales > 1/((T A) 
this transformation leaves only random backward scat- 
tering term of strength yitg + u>q/(4cto) in the action, 
which results for sufficiently small K in the Al-phase. 
Thus, to conlude we only find the three phases LL, MI 
and AI. 

Finally a word on the variational approach and replica 



symmetry breaking (RSB) which has been used in 11 1 
where the Mott glass phase has been found. In the vari- 
ational approach the full Hamiltonian is replaced by a 
harmonic one which leads to the decoupling of Fourier 
components (p q with different wave-vector q. Without 
RSB one obtains from this approach the results of per- 
turbation theory which is valid only on small scales. RSB 
gives the possibility of a further reduction of the free en- 
ergy. The results obtained in this way are exact only in 
cases when the coupling between different Fourier com- 
ponents is irrelevant and the physics is dominated by the 
largest length scale. Thus RSB is not an intrinsic prop- 
erty of the true solution of the problem, but a property 
of the variational approach. A good example is the re- 
lated problem of the interface roughening transition in 
a random potential 25]. The variational approach with 



RSB gives three phases: a fiat, a rough and a glassy flat 
phase [26j ]. where the rough phase has Flory like expo- 
nents. The functional RG (which takes the coupling of 
different Fourier modes into account) gives however only 
two phases: the flat and the rough phase, but at the tran- 
sition a logarithmically diverging interface width [25| . A 
similar situation may exist also in the present case al- 
though we are not able to calculate the properties of the 
AI to MI transition. 

To conclude we have shown for a one-dimensional dis- 
ordered Mott insulator (QJ, tracing back both the com- 
pressibility and the ac-conductivity to the kink energy 
Y, x of the electronic displacement field that an incom- 
pressible system has also a vanishing optical conductiv- 



ity. Thus we exclude the possibility of the existence of a 
Mott-glass phase. 
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